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Abstract. Chemical abundances of about fifteen ele- 
ments from oxygen to europium are measured in seven 
red giants of the two open clusters NGC 2360 and NGC 
2447. The effective temperatures of the giants are de- 
termined spectroscopically by taking advantage of their 
known masses (~ 2 M© in NGC 2360 and ~ 3 M© in 
NGC 2447) and bolometric magnitudes. 

The average iron abundances we obtain for the two 
clusters are [Fe/H]=0.07 for NGC 2360 and [Fe/H]=0.03 
for NGC 2447. 

Evolutionary stellar model calculations are performed 
in the mass range 1 — 4 M© in order to analyze the surface 
Na and O abundances predicted after the first dredge-up. 
The sodium abundance shows a well defined correlation 
with stellar mass in the 2 — 3 M© range. The agreement 
between our Na abundance determinations in NGC 2360 
and our model predictions at 2M© is very good. In con- 
trast, the overabundance in one of the three stars in NGC 
2447 exceeds that predicted at 3 M© by ~ 0.08 dex, which 
is significant compared to the observational error bars. 
The effects of core overshooting, convection prescription, 
metallicity and nuclear reaction rates on the Na surface 
predictions of our models are investigated. 

An oxygen deficiency relative to iron by 0.2 dex is mea- 
sured in our stars, in disagreement with our model predic- 
tions. Assuming that the Sun is 0.1-0.3 dex enriched in 
oxygen relative to neighbor stars could explain the dis- 
crepancy. 

Key words: Stars: abundances - Stars: fundamental pa- 
rameters - Stars: evolution - Stars: late-type - Open clus- 
ters and associations: individual: NGC 2360 - Open clus- 
ters and associations: individual: NGC 2447 



1. Introduction 

The comparison of chemical abundances observed at the 
surface of stars with those predicted by stellar model cal- 
culations offers a valuable tool to test our knowledge of 
stellar evolution. Among the most studied objects are red 
giants. Those objects are formed after the main sequence 
(MS) phase when hydrogen has been exhausted in the 
core. The core contracts, the outer layers expand, and the 
convective envelope deepens into the core. The consequent 
mixing of the ashes of hydrogen burning from the deep 
layers to the surface is called the first dredge-up process 
(1DUP). 

In most cases, the surface abundances after 1DUP are 
predicted to be sensitive to stellar mass. The observation 
of red giants belonging to given clusters presents in this re- 
spect a unique advantage since the age of the cluster, and 
thus of its individual stars, is known by comparing the 
cluster's distribution in the Hertzsprung-Russell diagram 
with predicted isochrones. The masses of the individual 
stars are then determined, enabling a thorough compari- 
son of their surface abundances with predictions. 

In this paper, we determine the abundances of oxygen, 
sodium and heavier elements in four red giants of the open 
clusters NGC 2360 and in three red giants of NGC 2447. 
All these giants belong to the clumps of their respective 
clusters. The initial masses of our stars are estimated to 2 
and 3 M© for those in NGC 2360 and NGC 2447, respec- 
tively. These masses turn out to be in a critical range for 
surface sodium enhancement after the 1DUP. Oxygen, on 
the other hand, is an interesting element to study since it 
has already been observed to be slightly deficient relative 
to iron not only in supergiants but also in giants of globu- 
lar clusters (Brown & Wallerstein 1991) and open clusters 



* Based on observations collected at the European Southern 
Observatory, La Silla, Chile (period 56) 

Table 3 is available only in electronic form from the Stras- 
bourg ftp server at 130.79.128.5 



(Luck 1994). Finally, the abundances of the iron-group el- 
ements are used to directly derive the metallicity of the 
two clusters. This direct technique, which contrasts with 
the classical, but less reliable, technique based on color 
indices and photometric calibrations, has so far been used 
only in a sm all number of open clusters (about 20, see 
Strobel 1991 , Twarog et al. 1997 ). Our observations en- 
able to enlarge that list. 
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The observational material is presented in Sect, 
the abundance determination procedure in Sect, 
metallicity is derived in Sect. ^, and the sodium and oxy- 
gen abundances analyzed in Sects. |5| and ^], respectively. 
Conclusions are drawn in Sect. ^. 



2. Observational Material 

Table [l] summarizes some parameters characterizing the 
two open clusters NGC 2360 and NGC 2447 studied in this 
paper. The distance modulus has been re-determined us- 
ing Mermilliod's database ( 1995 , 1999 ) and the age is de- 
termined with the theoretical isochrones of Schaller et al. 



( |1992|) . They lead to turn-off masses of 1.98 and 2.75 M 
for NGC 2360 and NGC 2447, respectively. 



Table 1. Adopted parameters for NGC 2360 and 2447, in- 
cluding the iron abundance [Fe/H] and metallicity [M/H] 
obtained in this paper (Sect. ||). The mass M tU rn-off at the 
turn- off (defined as point T in Fig. 1 of Maeder & Meynet 



199l| ) is given too. 





NGC 2360 


NGC 2447 


a (1950) 


7'T5.5 m 


7 h 42.5 m 


5 (1950) 


-15°32' 


-15°32' 


m v — My 


10.40 


10.25 


E(B - V) 


0.07 


0.04 


Ay 


0.22 


0.13 


Distance [kpc] 


1.086 


1.057 


Log Age 


9.06 


8.65 


Mturn-off [M ] 


1.98 


2.75 


[Fe/H] 


0.07 


0.03 


[M/R] 


0.10 


0.05 



The seven red giants studied in this paper have been 
selected on the basis of their constant radial velocity, as 
determined by CORAVEL observations, in order to ensure 
that they are single stars. This point is important since the 
light of a secondary component could bias the equivalent 
widths used to estimate the abundances of the primary. 

For each star, one echelle spectrum has been obtained 
by DR in November 1995 using the EMMI spectrograph 
attached to the NTT 3.5m telescope at the European 
Southern Observatory in La Silla, Chile. The grism #5 
and the echelle grating #10 were used, yielding spectra in 
the wavelength range from 4050A to 6650A with a resolv- 
ing power of R = 28 000. The detector was the ESO CCD 
#36, with 2048 x 2048 pixels. The signal-to-noise ratio 
per pixel of the extracted spectra varies from 54 to 226, 
depending on the wavelength range and exposure time. 

The journal of the observations is given in Table 0, 
together with the main characteristics of each stars (iden- 
tification, V magnitude, Geneva [B — V] color index, S/N 
ratio, Julian date, exposure time, radial velocity derived 



from our spectra and average RV from CORAVEL ob- 
servations) . The radial velocities are determined from our 
spectra by cross-correlation with a CORAVEL-type bi- 
nary mask optimized for the K0 spectral type. They are 
compatible with those determined earlier from CORAVEL 



observations by Mermilliod & Mayor (1989, 199C), con- 
firming that the red giants of our sample are very probably 
single. 

Observations of two IAU RV standards, HD 66141 and 
HD 80170 (of which Udry et al. |1999| give the precise RV 
value from CORAVEL observations), show that our re- 
sults can be trusted to within ±1.0 kms -1 . To reach this 
accuracy, we had to monitor the instrumental drift by tak- 
ing a few short Th-Ar calibration exposures during the 
night and correlate them with a longer calibration expo- 
sure taken in the preceding afternoon; then we linearly in- 
terpolated the instrumental velocity shifts for the epochs 
of the science exposures. 



3. Abundance determination 

For the reduction of the spectra, we used the Inter-Tacos 
software of Geneva Observatory]^]. The spectra were di- 
vided by an average of six flat-field exposures, after the 
background had been fitted by 2-d polynomials and sub- 
tracted. Then, they were extracted using a simple addition 
of the intensities along a virtual slit of 15 pixels. The wave- 
length calibration was performed with the same software, 
using Th-Ar spectra taken in the afternoon preceding the 
night. The rms scatter of the residuals around the fitted 
dispersion relation was 3.1 mA, which is quite good for this 
resolution. After rcbinning to a constant wavelength step 
of 0.03 A, which is over-sampled, we filtered all spectra 
by eliminating the high frequencies of their Fourier Trans- 
form. For that purpose, we used a FFT with a window of 
65% of AA (AA = A/R). The 55 orders were individually 
normalized to the continuum, using an interactive Super- 
mongo procedure which interpolates between the contin- 
uum points by a 3rd degree spline. 

3.1. Abundance determination 

The list of lines and oscillator strengths used for our abun- 
dance determinations is the same as that adopted by Bo- 
yarchuk et al. ( |1996[ ). The photospheric lines are chosen 
in wavelength zones avoiding the telluric absorption lines. 
The equivalent widths below 10 mA are rejected because 
they are too sensitive to noise and normalization. Table 3 
(available only in electronic form) summarizes the line pa- 
rameters and measured equivalent widths. For Mgi, not 
included in the list o f Boyarchuk et al., we adopt the log gf 
values of Thevenin fll989j |l990[). 



1 The Inter-Tacos software is the same as that used for the 
Elodie and Coralie spectrographs. 
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Table 2. Observational data of the 7 red giants and signal-to-noise ratio per pixel of their extracted spectra. The V 
magnitude and [B — V] index are those of Geneva photometry. The adopted numbering is that of Becker et al. (1976), 
but we also give the numbers in the system of Eggen (1968) for NGC 2360. There is no other numbering system for 
NGC 2447, but the map of this cluster is available on the web (Mermilliod 1999). The radial velocities of the 7 red 
giants and of two reference stars are also given, as determined in this work ("EMMI") at the indicated Julian dates. 
The exposure times are also listed. The CORAVEL RV values ("COR.") of the cluster stars are average ones from 
Mermilliod & Mayor (1989, 1990), while those of the reference stars (IAU standards) are taken from Udry et al. (1999). 
The errors quoted for the CORAVEL observations are related to the average RV values; the last column ("N") gives 
the number of CORAVEL measurements. 



Star 


DM or 


V 


[B-V] 


S/N ratio around 


JD 


t-exp 


V r [k 


ms- 1 ] 


N 




Eggen No 






5 150 A 


6600A 


-2450000 




EMMI 


COR. 


COR. 


2360-7 


8 


11.087 


0.294 


54 


89 


30.741 


1800 


28.53 


27.2±0.2 


5 


2360-50 


67 


11.082 


0.311 


71 


108 


30.767 


2700 


27.65 


27.2±0.3 


4 


2360-62 


81 


11.272 


0.251 


73 


146 


31.730 


3600 


28.05 


27.3±0.5 


4 


2360-86 


110 


10.787 


0.309 


114 


181 


31.822 


3600 


26.72 


27.3±0.3 


4 


2447-28 


-23° 6102 


9.849 


0.226 


112 


170 


30.847 


1800 


19.98 


21.2±0.2 


3 


2447-34 




10.123 


0.197 


143 


236 


32.755 


3000 


23.30 


22.1±0.2 


3 


2447-41 




10.031 


0.204 


124 


226 


32.840 


2700 


21.40 


21.5±0.2 


3 


HD 66141 


HR 3145 


4.406 


0.637 






30.833 


90 


71.43 


71.6±0.3 
















33.827 


90 


73.10 






HD 80170 


HR 3694 


5.310 


0.510 






31.799 


180 


1.38 


0.5±0.2 





The abundances are determined using the MOOG code 
(C. Sneden, Texas University) in the "abfind" mode, com- 
bined with atmospheric models of Kurucz ( |1996 ). The 
Kurucz models are interpolated by cubic splines in the 
T e ff and logg plane (solar abundances are assumed as 
far as the atmospheric structure is concerned). For some 
lines for which we have reasons to suspect that blends 
would make the equivalent width method risky (oxygen 
and europium), we compute a synthetic spectrum includ- 
ing a large number of neighboring lines, and determine the 
abundance of the considered element by best fitting the 
synthetic spectra to the observed one (the 'best' fit being 
evaluated by eye estimate). In that case, the MOOG code 
is used in its "synth" mode and the lines are taken from 
the VALD database (Piskunov et al. |l995| ). 

The measurement of abundance deviations as small as 
0.15 dex requires good spectra with signal-to- noise (S/N) 
ratios above 100. Such a S/N is reached at almost all wave- 
lengths in all our spectra, except near the [0 1] line in three 
stars in NGC 2360 (star numbers 7, 50 and 62). No O 
abundance is presented for those stars. 

3.2. Atmospheric parameters 

The determination of the atmospheric parameters logg 
(gravity), T c ff (effective temperature) and Vt (microturbu- 
lent velocity) is crucial for a good abundance estimation. 
The fact that the distance to the open clusters and the 
masses M of the individual stars are known provides a 
relation between log g and T e f j given by 

log 5 = -12.51 + logAf + 41ogT cfr + 0.4(M V + BC), (1) 



where M is expressed in solar mass, My is the absolute 
visual magnitude in V and BC is the bolometric correc- 
tion. We assume Afboi© = 4.75, and take the bolometric 
corrections from Flower (1977) (which are very close to 
the theoretically derived ones of Bessell et al. 



1998) 



The atmospheric parameters T e s and Vt of each star 
are estimated from an iterative procedure based on the 
method used by Boyarchuk et al. ( 1996 ), but modified to 
take into account the known bolometric magnitudes of the 
stars. 

Let us consider a given star. Starting with an initial 
guess of Vt — 1.6 kins" 1 , we construct a diagram (Fig. |IJ) 
plotting the abundances of the iron-peak elements com- 
puted by MOOG from our observed spectrum, as a func- 
tion of the effective temperature (log g is computed for 
each T e g value through Eq. 1). If all iron-peak elements 
had a solar abundance distribution, then there should be 
a value of T c g at which their abundances (normalized to 
solar) are identical. In practice, we choose (by visual esti- 
mate) the value of T c g at which the spread of abundances 
is minimal. Then, for this T e g- value, the usual diagram 
showing the Fe I abundances versus equivalent widths is 
used to adjust the best velocity of microturbulcnce V,. 
Finally, the whole procedure is reiterated in order to get 
mutually consistent parameters. 

The value of Vt found after the convergence of the iter- 
ative procedure differs by no more than 0.15 kms -1 with 
the first guess, which is consistent with the formal stan- 
dard deviation of 0.11 kms -1 for that parameter, obtained 
from the uncertainty on the regression line of Fe 1 abun- 
dances versus equivalent widths. 

This iterative procedure is applied to each of our stars. 
The resulting atmospheric parameters are summarized in 
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Fig. 1. Abundances of the iron-group elements as a func- 
tion of effective temperature for the star NGC 2447-41. 
The value of log g is given by Eq. [| at each effective tem- 
perature. The vertical line shows the adopted T e s value. 



Table [I| The derived effective temperatures agree to better 
than 100 K with the classical one based on the lack of 
correlation between abundance and excitation potential, 
the agreement improving with S/N ratio. 

Finally, the sensitivity of our abundance determina- 
tions to the adopted atmosphere models is tested by re- 
computing some abundances with the models of Bell et 
al. (1976). They lead to similar results, in the sense that 



abundances given by individual lines are the same within 
±0.07 dex after allowance for a slight systematic offset: 
the models of Bell et al. systematically lead to a 0.04 dex 
underabundancc compared to Kurucz's models (for lines 
with an equivalent width W\ < 100 mA; the offset is larger 
for stronger lines). In view of the relatively good agree- 
ment between both sets of models, we chose the Kurucz 
ones because they are more recent and easier to use. 

The surface abundance of a given element is deter- 
mined by adjusting the equivalent width(s) of its line(s) 
predicted by synthetic spectra to those measured in the 
observed spectra. Actually the MOOG program computes 
an abundance for each line of a given element, and derives 
a mean abundance value by averaging the abundances as- 
sociated to each of those lines. In computing the synthetic 
spectra, the formation of a dozen molecules (among them 
C 2 , CH, CN, CO, MgH, TiO) is considered. 

In the case of Eu, the only usable line is blended. We 
therefore determine its abundance by a visual match on 
a synthetic spectrum which includes all possible blends. 
This element is marked with an asterisk in Table 3. The 
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Fig. 2. Microturbulent velocity as a function of surface 
gravity. Full dots: NGC 2360; full triangles: NGC 2447; 
open dots: field giants of Boyarchuk et al. (1996); open 
triangles: cluster giants and supergiants of Luck (1994) 
with 4500 < T cS < 5500 K. Full line: regression line fit- 
ted to our seven stars assuming the same error on both 
axes; dotted line: regression line assuming no error on the 
abscissa. 



same technique is used for oxygen, discussed in Sect. [| 
The A6645.1 Euli line does not seem to be significantly 
affected by hyperfine structure, since its equivalent width 
in the Sun provides an abundance perfectly consistent 
with the meteoritic one, using precise experimental os- 
cillator strengths (Biemont et al. 1982). Our log gf value 
adopted for this line is only 0.03 dex below the one used by 
Biemont et al., so it would lead to a solar Eu abundance of 
0.53 instead of 0.50 dex, while the meteoritic abundance 



is 0.54 ± 0.01 dex (Grevesse & Noels |1993[) . 

The elemental abundances found in this work are listed 
in Table f§ 

3.3. Vt versus \ogg 

For information, we plot in Fig. |^ (filled circles and trian- 
gles) the microturbulence velocity deduced for our stars 
as a function of the surface gravity, together with the 
measurements of Boyarchuk et al. (1996) for five G7-K0 
field giants (open circles, their sixth star is not consid- 
ered because it is much cooler, of M0 type) and those 
of Luck (1994) for clusters' supergiants and giants with 
4500 < T off < 5500 K (open triangles). The combined 
data reveals a positive correlation between V4 and logg. 
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Considering only our seven points, the usual corre- 
lation coefficient p = —0.697 gives a t-test value of 
only —2.17, implying a significance level of about 90 
percent, while the Spearman rank correlation coefficient 
PSpearman = —0.786 gives a slightly better t-test of —2.84 
and a significance level of 96 percent. For the usual LSQ 
fit assuming no error on the abscissa and a constant error 
on the ordinate, we find 



V t = (-0.63 ± 0.29) log 5 + (3.37 ± 0.80) 



(2) 



with an rms scatter of the residuals of 0.083 kms" 1 . For a 
fit where errors are assumed to be the same on both axes, 
the relation would be 



V t = -0.861og 5 + 4.01 



(3) 



with an rms scatter of the residuals of 0.088 kms -1 . The 
scatter of the residuals is entirely compatible with the in- 
ternal standard deviation of Vt give n abo ve. It is smaller 
than that of t he da ta from Luck ( 1994 ) and from Bo- 
yarchuk et al. ( 1996 ), partly due to the better homogene- 
ity in metallicity and effective temperature. 

4. Metallicity 

It is difficult to derive an error estimate on the abun- 
dances. If we assume that each star has a solar distribu- 
tion of heavy elements, then an error estimate is given by 
the dispersion of the heavy element abundances around 
the mean abundance (relative to solar). The abundances 
of all measured elements are shown in Fig. || for each of 
our stars, with the mean iron-peak element abundance 
indicated by a dashed line in each panel. From visual in- 
spection (and excluding Ni which seems underabundant in 
some of our stars, probably only because of some system- 
atic error in the oscillator strengths), we estimate the error 
on the abundance of a given element to be ±0.1 dex. Ta- 
ble U summarizes the abundances derived in all our stars 
together with the number of lines used for their determi- 
nation. 

The metallicity of each star is defined as a weighted 
average of the abundances of the iron-group elements, 
the weight being the number of lines. Since these weights 
are roughly the same for all stars, each cluster's metallic- 
ity is defined as the arithmetic mean of the metallicities 
of its stars. The values found are [M/H]=0.10 for NGC 
2360 and [M/H]=0.05 for NGC 2447. The iron abun- 
dances are slightly smaller: [Fe/H]=0.07 for NGC 2360 
and [M/H]=0.03 for NGC 2447. Those are reported in 
Table g 

5. Sodium 

5.1. Abundances 

The sodium overabundances [Na/Fe] derived from our 
spectra are summarized in Table for our seven stars. 




W x [mi] 

Fig. 4. NLTE eff ects o n the Nai A6154.22 line according 
to Gratton et al. (1999). Full dots linked with continuous 
curves are for models with T c g = 5000 K and the three 
surface gravities logg = 1.5,3.0,4.5. Open squares linked 
with dotted lines: T e g = 6000 K and same log g values. The 
metallicity is solar ([M/H]=0.0) in all cases. The shaded 
area defines the range of equivalent widths observed in our 
sample of seven giants. 



They are small in NGC 2360, the oldest of the two clus- 
ters, with [Na/Fe]= 0.07 - 0.13 dex, and amount to 0.15 
- 0.22 dex in NGC 2447. 

The error in the Na overabundance [Na/Fe] due to 
the uncertainties on T e gf, log<? and Vt are not large be- 
cause Na I and Fe I have similar behavior s with respect to 
those parameters (e.g. Table 9 of Luck 1994 ) such that 
their effect on the [Na/Fe] ratio cancels out. The errors 
on the [Na/Fe] ratio due to the continuum placement, 
on the one hand, and to the photon noise, on the other 
hand, are estimated in the following way. For each star, 
the equivalent widths of the two Na I lines are measured 
five times, each time after having renormalized the spec- 
trum of the relevant region. The standard deviations of 
the equivalent widths are then computed from these data 
and propagated into standard deviations of abundances 
using the MOOG code. The resulting dispersion on the 
Na abundance amounts to 0.02-0.05 dex, depending on 
the S/N ratio. A synthetic spectrum of a typical red gi- 
ant (with T cS = 5200 K, log g = 2.77, [M/H]=0.0 but 
[Na/Fe]=+0.15) is then produced, and a Gaussian noise 
added to it for four representative S/N values, five spectra 
being produced independently with the same S/N ratio. 
The equivalent widths of both Na lines are measured (by 
a Gaussian fit) for these 20 spectra and their standard 
deviation computed for each S/N ratio. The deviations 
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Fig. 3. Abundances relative to solar for the 7 red giants. The dotted line represents the solar abundances, the dashed 
line is fM/Hl. defined as the a.vera.ee abundance of the iron-nea.k elements weighted bv the number of lines. 
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are found to lie between 2 and 3.3 percent for S/N ratios 
between 230 and 90. They are translated into abundance 
errors (after dividing them by v2 since there are two lines) 
and added quadratically to the errors due to the contin- 
uum position. Finally, we assume that the error on the 
continuum position for Fe I is similar to that for Na I and, 
admitting that they are independent of each other, add 
them quadratically to the total error on the Na abun- 
dance^] (in fact there is probably a correlation between 
the continuum placement for the Fe lines and that for the 
Na lines, but neglecting it only results in an overestimate 
of the error, so that we stay on the safe side). The result- 
ing estimated errors amount to about 0.03-0.07 dex, the 
lowest ones pertaining to NGC 2447. 

Among the systematic errors which may affect the Na 
and Fe abundances, some are negligible for the [Na/Fe] 
abundance ratio because of their mutual cancelation, as 
mentioned above. This is not the case, however, for the 
oscillator strengths adopted in the synthetic spectra. If the 
log gf values of the two Nal lines are slightly in error while 
those of Fel are statistically correct, for instance, then the 
absolute values of [Na/Fe] would be wrong, though their 
relative values (i.e. their differences) would remain valid. 

Finally, let us consider the systematic errors due to the 
assumption of local thermodynamic equilibrium (LTE) in 
the MOOG program. No n-LTE (NLTE) calculations per- 
formed by Gratton et al. (199E) in atmospheric conditions 
relevant to red giants show that the NLTE effect on the 
Nal A6154.22 line (one of the two lines used in this pa- 
per) strongly depends on surface gravity and slightly on 
effective temperature. This is illustrated in Fig. |4], which 
reveals that the NLTE correction should be small in the 
2.5 < logg < 3 range characterizing our stars. Indeed, the 
average NLTE corrections are found to amount to 0.006- 
0.035 dex for the stars in NGC 2360 and of 0.032-0.048 dex 
for NGC 2447 using Table 11 of Gratton et al. (quadrati- 
cally interpolated in log g and linearly in T c g, and assum- 
ing similar corrections for both A6154.22 and A6160.75 
lines of Na i) . A similar calculation for iron^j leads to NLTE 
corrections for Fe of 0.015-0.021 dex for NGC 2360 and 
0.027-0.028 dex for NGC 2447. The resulting NLTE effects 
on [Na/Fe] thus ranges between -0.013 and 0.014 dex for 
NGC 2360 and between 0.005 and 0.020 for NGC 2447. Of 
course, these values are only approximate since they are 
derived from only one iron line assumed to be representa- 
tive of the 55 lines observed for that element, but they do 
suggest that the errors brought by our LTE approxima- 
tion are much smaller than the above mentioned random 
error bars. 



2 the error on the Fe abundance due to photon noise is ne- 
glected because of the large number of lines of this element 

3 We take the theoretical Fel high ex citatio n line (EP = 
4.1 eV) with W\ = 60 mA (Gratton et al. 1999) as representa- 



Surface S3 Na overabundance 
after first dredge-up 

NGC 2447 




M (M e ) 



Fig. 5. Na overabundance versus stellar mass. The rect- 
angles represent the abundance determined in our seven 
red giants with the estimated error bars. The theoretical 
curves are obtained at metallicity [Z=0.05]. The continu- 
ous line joining the full dots is obtained with standard re- 
action rates and without overshooting, while that joining 
the open dots takes overshooting into account. The dotted 
lines show the predictions of a standard model (without 
overshooting) with extreme reaction rates still compatible 
with the uncertainties. 



As a conclusion, the main abundance errors, be- 
sides the possible systematic errors due to the oscillator 
strengths^ , are due to the equivalent width and continuum 
measurements, which amount up to 0.07 dex. The error 
bars on the masses of each star, on the other hand, are 
estimated by considering a 0.05 dex error on the age logi 
of the clusters. The surface sodium overabundances as a 
function of stellar mass are shown in Fig. || by rectangles 
taking into account the above mentioned uncertainties. 

5.2. Predictions 

Sodium production during H-burning results from the 
transformation of 22 Ne into 23 Na by proton capture. This 
reaction occurs very efficiently at the temperatures charac- 
terizing the core of MS stars (see Appendix A of Mowlavi 
1999). First dredge-up then mixes some of the synthe- 
sized Na from the deep layers to the surface. This sce- 
nario is confirmed by the observation of sodium overabun- 
danc es at t he surface of ma ny gia nts and supergiants (e.g. 
Luck fL994l Boyarchuk et al. I1996L Takeda & Takada-Hidai 



1994) 



The surface Na overabundance predicted by stellar 
model calculations (without core overshooting) as a func- 



The log gf values of the Na I lines given by Boyarchuk et al. 



tive of our data. Taking the low excitation line (EP = 0.0 eV) 
or W\ — 90 mA would not change the results for the NLTE 
corrections. 



(|1996j) are 0.05 dex higher than those of T heven in ( |l99C| ), but 



0.04 dex lower than those of Wiese et al. ( |l969j ); this suggests 
that the systematic error on the log gf values - hence on the 
abundances - cannot be much larger than 0.05 dex. This is 
consistent with the error estimate of Thevenin. 



8 S. Hamdani, P. North, N. Mowlavi, D. Raboud, J.-C. Mermilliod: Chemical Abundances in 7 Red Giants 



0.6 



0.4 



+ 4- 



•0.2 



IB 

y 



+ XGC 2447 

% NGC 2360 



-0.2 



log(g) [cgs] 



Fig. 6. Na overabundance versus log g. The crosses repre- 
sent yellow supergiants and giants in clusters measured by 
Luck 1994), the open triangles represent red giants mea- 
sured by Boyarchuk (1996) and the full symbols repre- 
sent the seven red giants of this study (full dots: NGC 
2360, full diamonds: NGC 2447). The two crosses at high 
log<? are dwarfs, while the two supergiants with a nega- 
tive [Na/Fc] value arc among the coolest in Luck's sample 
(T off = 4000 K). The short solid line on the right of our 
observations corresponds to predictions of clump stars of 
2.2 to 2.8 M Q computed with core overshooting. 



tion of stellar mass is shown in Fig. || by filled circles 
connected with solid line. The models have a metal- 
licity 0.05 dex above solar (which is the metallicity of 
NGC 2447) and are followed from the pre-MS up to the 
completion of the 1DUP . The stellar evolution code is 
the same as in Mowlavi ( 1999 ), exce pt that the NACRE 
reaction rates (Arnould et al. 1999] ) are used to follow 
the nucleosynthesis and that the formalism of Canuto, 
Goldman & Mazzitclli (1996) is used to describe the en- 
ergy transport in convective zones. The surface sodium 
abundance is seen to be very sensitive to stellar mass in 
the 1.5 — 3M Q range. It increases from no Na enhance- 
ment below ~ 1.5 M© to an overabundance of 0.15 dex 
at ~ 3 M Q , and keeps this value for stellar masses up to 
6M0 before increasing again with stellar mass (see, e.g., 
Mowlavi 



1998) 



The sensitivity of those predictions to core overshoot- 
ing, convection prescription and stellar metallicity is ex- 
plored by computing extra models from the pre-MS to the 
1DUP. None of those parameters, however, turns out to 
have a significant impact on the surface Na abundance. 
Models with core overshooting (with an extra-mixing ex- 
tent of 0.20 times the pressure scale at the core boundary) 



predict a 0.02 dex enhancement (open circles connected 
with solid line in Fig. ^J) compared to predictions without 
core overshooting. Increasing the metallicity by 0.17 dex 
does not change the surface [Na/Fe] prediction after 1DUP 
by more than 0.01 dex. And using the mixing length the- 
ory (with a mixing length of 1.5 times the pressure scale 
height) instead of the Canuto, Goldmann & Mazzitelli for- 
malism does not change the surface abundance predictions 
within 0.001 dex. 

Let us now explore the uncertainties linked to nuclear 
reaction rates. Both the 22 Ne (p , 7) 23 Na and p-capture 
reactions on 23 Na ar e still subject to large uncertainties 
(Arnould et al. 1999 ). In order to assess their impact on 
our surface Na abundance predictions, several models are 
recomputed from the pre-MS up to the 1DUP with the 
upper/lower limits for the rates provided by the NACRE 
compilation (cf. Arnould et al.), as appropriate to mini- 
mize/maximize 23 Na production. The results in the 'min- 
imal' and 'maximal' 23 Na production cases are shown in 
filled circles connected with dotted lines in Fig. ||. They 
reveal a variation in the surface Na abundance predictions 
of up to 0.08 dex relative to the 'nominal' case where the 
recommended NACRE rates are used. Nuclear reaction 
rate uncertainties thus dominate the uncertainties associ- 
ated with stellar metallicity and convection prescriptions 
for sodium predictions in red giants. 

5.3. Discussion 

Figure || shows a very good agreement between our Na 
abundance predictions in 2 M Q model stars and those ob- 
served in NGC 2360. The nominal predictions of the 3 M Q 
model star, on the other hand, seem a little too low com- 
pared to the abundances observed in NGC 2447. The pre- 
dictions in the maximal case of 23 Na production would 
fit the highest Na abundance measured among the three 
stars observed in NGC 2447. Those predictions, however, 
would not be compatible with the Na abundances mea- 
sured in the 2 M red giants of NGC 2360. Figure | thus 
suggests that the nominal Ne-Na reaction rates should 
not be too much altered, if at all. The solution to the 
discrepancy between our observed Na abundance obser- 
vations in NGC2447 and predictions should be found in 
other(s) mechanism(s) such as, possibly, meridional mix- 
ing induced by stellar rotation. Further theoretical and 
observational investigations should be performed before 
being able to draw a firm conclusion. 

Finally, let us mention that the positive sodium abun- 
dance - stellar mass correlation translates, at a given ef- 
fective temperature into a sodium abundance - surface 
gravity anti-correlation (or sodium abundance - luminos- 
ity correlation). This is well known in the literature, and 
shown in Fig. || where our data are displayed together 
with those of Luck (1994) and Boyarchuck et al. (1996). 
The dependence on effective temperature is small (Luck 
|l994|) . The computation of our 2.2 to 2.8 M models with 
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core overshooting is carried on up to the clump in the core 
helium burning phase. The clump is defined as the point 
where the stellar luminosity (gravity) reaches its minimum 
(maximum) value after core helium ignition. The sodium 
abundance (which is not altered between the 1DUP and 
the clump) predictions for those clump models are shown 
by a solid line in Fig. 6. 



6. Oxygen 

6.1. Abundances 

The O abundance is determined from the 6300. 3lA line 
of Oi, which is available in our spectra with a S/N ratio 
above 100 only for the stars of NGC 2447 and for one 
of the stars in NGC 2360. This line is slightly blended 
with a line of Sen at 6300. 68A, but both lines are fairly 
well separated in our spectra (see Fig. ||). Fortunately, 
no significant telluric line is spoiling the O and Sc lines, 
though some are present 1 A away or more. Using the 
equivalent width of the 6300. 3lA line, the O abundance 
is found to be slightly deficient with respect to iron, with 
[0/Fe]= -0.1 to -0.3 dex (see Table @). 

The Oi line at 6300. 3lA, however, may be blended 
with another line at 6300. 33A from Nil. The logg/ of 
the Nil line is still poorly known. If we take the only 
data available in the recent literature, i.e. log gf = —1.737 
and EP 



4.266 eV proposed by Kurucz (1994), and 
if we assume [Ni/Fe]=0.00, then the oxygen abundance 
predicted by synthetic spectra reaches -0.6 dex on aver- 
age (cf. [O/Fe]* in Table 0). A smaller value of log gf 
is however suggested from the analysis of the solar spec- 
trum. Lambert ( 1978Q estimates the width W\ of the Nil 
A6300.33 A line in the solar spectrum to be below 0.1 
- 0.5 mA. The oscillator strength derived from the solar 
Ni abundance would be ten times smaller than Kurucz's 
even with W\ = 0.5 mA. The true log gf of this Nil line is 
thus very probably smaller than -2.74 dex. Therefore, the 
blend cannot affect the oxygen abundance determination 
by more than a few hundredths of dex. 



6.2. Predictions 

The surface O underabundance predictions after the first 
dredge-up are shown in Fig. [?] as a function of stellar mass. 
Models with (open circles) and without (filled circles) core 
overshooting during MS both predict very little variation 
of the surface O abundance in red giants (depletion by less 
than 0.04 dex). This results from the fact that only the 
deepest stellar layers are affected by the N-0 cycle of H- 
burning. As a result, the dilution factor of the envelope's 
oxygen is rather small when the envelope extends down to 
the inner regions. 



fc. -0.3 



Surface le 
underabundance 
after first dredge-up 



l 

0.9 

0.8 



M (M s ) 

Fig. 7. Same as Fig. 0, but for oxygen. 
6.3. Discussion 

The predicted surface O underabundances fall much below 
the observed ones (by 0.1 - 0.3 dex if Ni is not taken 
into account and by more than 0.5 dex with Ni). This 
discrepancy is unexplained so far. 

Such a deficiency of O in red giants has already been 
report ed in the literature (e.g. Luck 1994 ). Interestingly, 
Venn J199SD argues that the Sun presents an overabun- 
dance of oxygen of about 0.3 dex relative to the galactic 
A supergiants, to B stars in the solar neighborhood and 
in Orion, and to nebular abundances in the Orion nebula. 
The conclusion that the solar oxygen abundance is higher 
than what should be expected by galactic chemical evolu- 
tion is also reached, at least qualitatively, from the sample 



of G dwarfs of Edvardsson et al. (1993), which show an 
average oxygen abundance of [O/Fe]~-0.1 at [Fe/H]>-0.1. 
If this empirical fact is confirmed, it would explain in a 
simple way the oxygen underabundances found in our gi- 
ants, at least if we neglect the contribution of the Nil line 
in the synthetic spectra. 

7. Conclusions 

The abundances of about fifteen elements in seven red gi- 
ants of NGC 2360 and NGC 2447 are derived with a global 
uncertainty of about ±0.1 dex (essentially due to photon 
noise, normalization of the spectra and uncertainties re- 
lated to atomic line parameters). From those abundances, 
the following conclusions are drawn. 

The metallicity [M/H] of the two open clusters is di- 
rectly derived from their iron-group elements. It amounts 
to 0.10 for NGC 2360 and 0.05 for NGC 2447. 

The sodium abundance reveals a positive correlation 
with stellar mass in the 2 — 3 M Q range, as expected from 
model predictions. Quantitatively, the surface 23 Na abun- 
dances are in good agreement with those predicted after 
1DUP for 2M Q stars (i.e. for stars in NGC 2360), but 
not at 3M (NGC 2447). The disagreement in NGC 2447 
reaches 0.08 dex, which we consider marginally significant 
given our estimated uncertainties on the Na abundances. 
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NGC 2447-41 
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5200 K, logg-2.65, v,-1.70 km/s 
= -0.67, [So/H]=-0.31, [Ni/H] = 0.05 
FWHM(g. smooth) = 0.24 A 
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Fig. 8. Measured spectrum of the star NGC 2447-41 
(dots) , showing how the [O i] line used for the abundance 
determination is slightly blended with an Sen line. The 
continuous line is the synthetic spectrum computed with 
MOOG, including the [Oi] 6300, Sen and Nil lines. The 
abundances adopted for the fit are given on the plot. The 
measured spectrum has been slightly smoothed using the 
Fourier method, with a window of 0.15 A. 

The excess of Na observed in two of our stars of NGC 2447, 
if confirmed, cannot be explained by current stellar mod- 
els. Other physical mechanisms, such as stellar rotation, 
should be investigated in the future. More observations in 
open clusters of various ages (and thus turn-off masses) 
are also required to increase the small number statistic of 
this study. 

The oxygen abundance is very deficient (by 0.2 to 
0.5 dex depending on the true gf value of the Ni I A6300.34 
line) relative to iron in all our four red giants for which 
the oxygen abundance could be measured. This is not 
compatible with stellar model calculations which predict 
an oxygen deficiency by at most 0.03 dex after the first 
dredge-up. The high deficiencies can be partly explained 
if we assume that the Sun is slightly oxygen-rich relative 
to most stars in its neighborhood, as suggested by several 
authors. If not, then a new mechanism in stellar models 
should be thought of to account for this discrepancy. 
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Table 4. Fundamental parameters of the 7 red giants and mean abundances relative to the Sun. For each species, 
the second line gives the rms standard deviation of the abundances obtained from individual lines. N is the number 
of lines measured. The oxygen abundances and [O/Fe] are derived from the equivalent width neglecting the possible 
contribution of the Nil A6300.33 line, while [O/Fe]* is computed from a synthetic spectrum including the Nil A6300.33 
line (with the gf value from Kurucz 1994). [Fe/H] is an average weighted by the number of Fei and Fen lines. [M /H] 
is a similar average, but on all iron-peak elements from Ca to Ni. 
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